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ESTIMATING EVAPOTRANSPIRATION 


UNDER NONHOMOGENEOUS FIELD CONDITIONS ' 


By 

2 

Leslie H. Parniele and Earl L. Jacoby, Ji\ 


SUMMARY 


Detailed micrometeorological measurements were 
ade during selected days throughout 1969 by using 
weighing lysimcter on an adjacent watershed that 
been planted with corn in order to test the 
^plication of commonly used meteorological 
ethods of estimating evapotranspiration under con- 
itions of rolling terrain. The experimental site was 
dected on a watershed that had a slope of 1 3 percent 
It that consisted of rolling-surface topography 
mated in the U.S. Depaitment of Agriculture’s 
orth Appalachian Experimental Watershed, 
oshocton, Ohio. 

With average hourly data being used, the Bowen 
itio energy budget method proved superior to the 
ability-adjusted aerodynamic methods in com- 
nring hourly lysimeter flux rates. 

The Bowen method appears to be capable of 


accurately estimating actual evapotranspiration 
under natural field conditions where fetch and 
surface homogeneity is less than ideal. Hourly com- 
putations in the combination equation for potential 
evapotranspiration agreed with the lysimeter for bare 
soil and for corn 75 centimeters high. 

However, using the log profile wind function in the 
combination equation was unsatisfactory for corn at 
its full height of 240 centimeters unless the roughness 
length parameter (zo) was adjusted as a calibration 
term. An unusually low value for Zo of 1 .0 centimeter 
resulted in the best fit of the combination equation 
with lysimeter flux rates for corn higher than 75 
centimeters. Techniques were presented for making 
probable-error analysis of evapotranspiration 
estimates through meteorological methods. 


INTRODUCTION 


The evapoiation from baie soils and the 
^ipotranspi ration from vegetated surfaces have 
3cn studied in detail using meteorological methods 
ipable of providing accurate, short-period flux 
lies. Under ideal conditions, meteorological 
ethods have been highly successful in providing 
Hailed information on the vapor transfer processes. 


iConuibutcd by llic Northeast Walerslicd Research Center, 
.gMCUtluial Research Service. U.S. Depart mcnl of Agncullure, 
fnivei.sily Park, P«Im in cooperation with the Agricultural 
xperinienl Station, Pennsylvania State University, 

^Pormcrly agricultural engineer with the Northeast Watershed 


However, these methods have not been tested 
adequately for the nonhomogeneous sui faces, soils, 
and microclimatic conditions encountered in many 
natural situations apart from the ideal infinite fetch 
and short-growing, well-watered crops that 
characterized most past verification studies 
Experimental studies are needed to test the 
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application of nietcoiologkal tlicoiy against the 
realities of evapolianspiiation fiom the nonkical 
cioppcd siu faces that arc more commonly 
encountered in agricuUiua! practices. 

Thus, the available meteorological methods should 
be tested under the realities ot evapotianspiiation 
fiotri cropped and heterogeneous sin faces. The 
instrumentation requirements and assumptions basic 
to a given model can be verified by comparing com- 
puted results against an independent measure, pre- 
ferably an accur ate weighing lysimeter. The purpose 
of this pLiblrcation is to present such experimental 
r esults under less-than-ideal site conditions. 

Energy Balance Method 

One of the most promising meteorological 
methods available is an energy balance approach 
known as the Bowen ratio method (2).^ The Bowen 
ratio, defined as the ratioof sensible heat flux (H) to 
latent heat flux (E), can be derived by dividing the 
eddy transport equation for sensible heat flux by that 
for latent heat flux 

c PK|^ OT/az)p 

B = " = 

E LeK^Oe/9Z)p 

Assuming = K y, an approximation of B can be 
obtained by using a finite difference form; 

B = T 

Lg (Ae) Ae (2) 

where T and e are temperature and vapor 
pressure measured over the same height interval just 
above the vegetation or baie soil surface. 

Estimates of E using the Bowen ratio appioach are 
then calculated from. 


Rn + G 


where R and G are the measured net radiation and 
soil heat flux rates. Equation 3 has been used to 
measure evaporation from homogeneous surfaces, 
assuming that no horizontal divergence of fluxes 
occurs between the levels of measurement and the 
evaporating surface. 


Midhc nurnhets in parentheses icfci to I Ucialiue Cited, p. 19, 


The advantages and disadvantages of the Bowen 
method arc discussed mother souiccs (5, 26), I anner 
(27) repoited close agreement between calculated 
evapoti’anspiration lates and lysiinctei “measured 
values in a humid climate with a wclhwatei ed crop In 
warm, arid le^ons, the Bowen method has been used 
by PTitschen (5, 6) successfully to determine tlie 
evapoiative flux for* short periods. These deter- 
minations can also be made for longer per iod values. 

Aerodynamic Approach 

Aerodynamic appi caches for estimating total 
water vapor flux from an evapor ating surface depend 
on assumptions that are well documented in the 
literature cited (2J, 26, 28, 22), Generally, these 
assumptions are: 

® Mean wind speed varies logarithmically 
with height. 

« Eddy shear stress is constant with height 

® Eddy diffusivity for water vapor (Kv)i^ 
equal to that for momentum (Km), 

The equation based on these assumptions most 
commonly used in evaporation measurement under 
adiabatic conditions is: 

n = .Lpk2(“l-‘-2) 

(In (Z2-D/Z,-D)|2 

(4) 

Equation 4 permits the prediction of latent heat flux 
rates from carefully measured profiles of vapor 
pr cssirie(c) and horizontal wind velocity (u). The dis- 
placement height (D) (D = + d) must be 

subtracted from the height value (Z) as mcasui'cd 
from ground surface, This adjustment is necessary 
because, as figure 1 shows, whei*c there is a tall,acro- 
dynamically rough vegetation, the logarithmic wind- 
profile distributions on which equation 4 is based ai e 
displaced upward from the ground surface. The 
parameter is a constant that is used to define the 
roughness characteristics of a given surface. It is 
determined by gr’aphically extrapolating the log- 
arithmic pi’ofilc of measured velocities downward to 
where wind speed theoretically should be zero, The 
parameter d is the height of the r’cfcrcncc plane 
displacement. 

The above approach is valuable only under ncai- 
adrabatic conditions. When surface heating or 
cooling produces a nonadiabatic temper atui c pr ofile, 
the straight-line reialronship of u versus I n (Z-D) no 
longer holds, especially under low-wind conditions. 
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hiijiij e 1 . —Schematic repi esenlation of wind velocity dblribiHion 
(U) with height above ground (Z) in a plant community with 
height (h) Wind velocity distribution is logarithmic with height 
above plants, and exponential with height within the plant 
community. Loganthmic piolilc extrapolation downward to 
/CIO wind .speed (u = o) is indicated by a dashed curve with an 
inteicepl of d Zq (/2) 


g (ax/az) 

Ri = 

(6a) 

where g is the constant of gravity, and where is 
the absolute tempeiature. Equation 6a can be 
approximated by 

Ri=(B(Z2-Z,)/T,b,)[(T2-T,)/{u2-Ui)2] 

(Ob) 

and is applicable to a mean height (z) between z ; and 
Z 2 , where AT and Au are measured. 

The coefficient b in equation 5 must be determined 
experimentally. Panofsky, Blackadar, and MeVehil 
(76) recommended a value of 18 after analysis of 
several data sources. With this b value, equation 5 is 
commonly called the “KEYPS” profile. 

Another empirical diabatic profile correction 
developed at Davis, Calif. {20) is; 

^ ^ {\ ±50 Ri)± 1/4, 

(7) 

where the -50 and -1 /4 apply in the unstable case, and 
the +50 and +1/4 apply in the stable case. 

Equation 4 can be modified by introducing the 
diabatic correction function \l/ , Thus, to correct 
evaporation flux calculations for all conditions of 
thermal stability, the following equation (26) is used 


Under conditions of appreciable temperature 
gradients because of surface warming or “free con- 
vection,” equation 4 will underestimate the flux rate 
since the added effect of buoyant turbulence has not 
been considered. Similarly, temperature inversions 
over cold surfaces damp out turbulent exchange, 
causing equation 4 to overestimate exchange rates. 

Several means have been proposed to adjust for tne 
influence of temperature gradients (23), One 
approach is to incorporate a form of diabatic 
correlation function that can then be related 
empirically to a stability parameter. 

The form of the relationship normally used is: 


where Ri is the gradient Richardson number, a com- 
monly used stability parameter. This non- 
dimensional number (Ri) is a quantitative ratio of 
thermal turbulence effects, compared with frictional 
turbulence effects that have the form of 


E = eLpk (Cj - 62) (U2 • 

[In (Z 2 - D) (Z, - D) [2 (,/,)2 (gj 


Equations for Potential Evapotranspiration 

Potential evapotranspiration is a micrometeoro- 
logical concept that is proposed to be the actual 
evapotranspiration under conditions of freely avail- 
able water and that is intended to provide a base line 
for computing actual evapotranspiration as a 
fraction of potential evapotranspiration. 

The combination method for estimating potential 
evaporation was developed by Penman from both 
energy balance and aerodynamic transport equations 
(77). Much attention has been given to the derivation 
and application of the combination method and its 
various versions in estimating potential evaporation. 
Many sources are available for more detailed infor- 
mation (7, 3. 7ft 77, 7ft 24, 27. 28, 31), 
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Van Bavel ( 3 !) described the development, 
limitations, and adequacy of the Penman com- 
bination concept in considerable detail. Moreiecent 
formulations of the combination equation are an 
improvement over the original Penman versions in 
that they contain fewer empirical constants or 
functions. Basically, the current accepted equation 
{ 31 ) is: 

s + 7 s Tg 

( 9 ) 

where e ^ and e a are the saturated and actual vapor 
pressure of the air at height Za above the given 
suiface In contrast with the original Penman 
formula: 

® Net radiation (Rj^ ) is actually measured, instead 
of relymgon an empirical estimate of this parameter. 

® Correction is madeforsoil heat nu.x(G) that had 
been previously oniited 

« Fy is derived from the log wind-profile theoiy of 
Thornthwaite and Holzman (i9)and referenced from 
the displacement height (D). 



« An empirical constant is not included as an 
adjustment for crop wetness or stomatal control. 

® Vapor pressure at the evaporating surface is 
equal to the saturation vapor pressure at surface 
temperature. 

The combination method for potential ET requires 
these assumptions: 

1. The eddy diffusivity for momentum equals that 

for heat and water vapor, = ICv 

2, The air temperature at point of measurement 
can be substituted without undue error for the 
average of the air and surface temperatures in 
evaluating 

r(s/s+7) 

By referring to equat ion 9, one canapply a method 
based on potential evapotranspiration to any surface 
with freely available water. The only differences in 
evaporation from a water surface, from bare wet soil, 
and from a well-watered vegetative surface would be 
due to variations in surface radiative properties and 


in aerodynamic transpoit characteristics of the layer 
through which water vapoi moves away from the 
evaporating surface, assuming that vapor pressure 
(Cs) is at saturation. 

In this form, no provision is made for differences in 
stomatal characteristics making one crop ‘Vetter” 
than another, even when well watered. Tests of com- 
bination equations 9 and 1 0 for open water, wet bare 
soil, and well-watered alfalfa (a “wet” crop) showed 
excellent agreement on both an hourly and on a daily 
basis under a variety of physical conditions in 
Phoenix, Ariz, (J/). 


Limitations of Meteorological Methods 

For routine field conditions, these meteorological 
methods could be impractical because of measuie- 
ment difficulties and inaccurate because o f the failure 
of the theoretical model to fit less-than-ideai con- 
ditions of a horizontally uniform, two-dimensional 
surface. Peterson ( 19 ) discussed the problems 
associated with flows encountered in heterogeneous 
terrain. His predictions indicate that the standaid 
techniques for estimating the eddy diffusivities oi 
exchange coefficients may lead to unacceptable 
errors in estimating the vertical fluxes of momentum, 
heat, and moisture over a heterogeneous surface. 
Many agricultural and natiiial surfaces simply can- 
not be regarded as homogeneous. 

Surface nonhomogeneity may not only affect the 
ratios of (K|^ /K„^) and (Kj, /K^,), but a Iso may create 
sampling problems in measui emenl of gradients and 
other terms with the required accuracy. 

When sensors are elevated to obtain better spatial 
sampling and maximum possible vertical differences, 
the profile measurements of wind, temperature, and 
vapor concentration are affected more and more by 
surfaces further upwind and less by the surface 
immediately below the sensor, At some point the 
measured profile gradients are not representative of 
processes that are occurring on the reference surface. 
Thus, a buffer area of the same type as the reference 
surface is required to reduce the possibility of 
divergence errors in measuring the vertical fluxes. 
This buffer area is referred to as "fetch” and is the 
minimum upwind distance to a significantly different 
surface from the point of measurement. 
Nonhomogenieiy influences fetch requirements in 
that the required height-to-fetch ratios should be 1 to 
50 for minor discontinuities and at least 1 to 1 00 or 1 
to 200 if the discontinuity is severe ( 26 ). These "ideal” 
fetch requirements are difficult if not impossible to 
satisfy under many field situations. 
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EXPERIMENTAL PROCEDURE 


The experimental site was developed on water** 
shed 109 at the U.S. Department of Agriculture, 
North Appalachian Experimental Watershed, 
Coshocton, Ohio. The watershed (fig. 2) lies just 
below the crest of a major ridge that extends through 
the area and that has an average slope of 1 3 percent m 
an eastward direction. Agricultural fields border the 
north, south, and west sides of the watershed; gravel 
roads travel se the north and east sides. The rolling 
terrain is far from an ideal homogeneous surface as 
required for the classical type of micro- 
meteorological studies. At the center of the watei- 
shed a micrometeorological sampling site for tern per- 
atuie, humidity, and wind profiles was established. 
Most energy budget data and further profile 
measurements of humidity, temperature, and wind 
were obtained at lysi meter Y102C, located 200 feet 
north of the cential wateished sampling site. The 
instrument trailei that housed the recording equip- 
ment was parked on the boundary of the field 
approximately 100 feet northwest of the lysiinctei 
site. 

The entire watershed, including lysimeter, was 
planted in corn on May 1, 1969, in 42-inch rows 
generally oi iented in a southwest-northwest direction 
parallel to the contour of the field. The prevailing 
winds in the area are from the southwest and provide 
the best fetch conditions foi the two sites. 

Complete micrometeorological measurements 
were made at the lysimeter site on May 2 and 22, 
June 26, July 24, August 27, and October Hand 16, 
1969, Micrometeorological measurements were not 
made at the midwatershed site during the first two 
test periods. The midwatershed profile data on 
June 26 were inaccurate because of instrumental 
difficulties, so that data at the midwatershed site were 
available only on July 24 and afterward. A summary 
of conditions for each data is given in table 1. 

The corn was harvested on October 1 5 and the soil 
disked to prepare for winter wheat planting. Thus, 
the data on October 16 were obtained on bare soil 
partly covered with litter from cornstalks. 

A complete description of the instrumentation and 
data-reduction techniques can be found in 
appendix B. A listing of the average hourly micro- 
meteorological data is listed in appendix C. 


Variation of D and Z 

The deterniinations of Dand that vveie used m 
equations 4 and 10 are two of the most demanding 
measurements in the profile methods, because each 
depends on accurate measurement of the second- 
order properties of the profile D and Zo must be 
detei mined under adiabatic situations because lack 
of stability in these paramenters may introduce errors 
{15), Since DandZo of ciopped surface may change 
several centimeteis with varying wind speed {12, 30), 
errors may also arise from assuming that D and 
are constant for a given surface. 

Miikammal (/5) has indicated that the 
determination of D is vei y sensitive to small errors in 
wind speed measurement He obseived that under a 
given set of conditions, an increase of only 1.0 
centimeter per second in the average speed of the 
middle level of a three-level wind system increased D 
by 50 percent, from 40 centimeters to 60. 

Error Analysis 

The meteorological estimates of E were basically 
indiiect measurements; that is, the unknown (E) was 
determined by using an equation that related it to 
some other measurable quantities such as radiation, 
temperature, wind speed, and vapor pressure. Each 
measured vaiiable had a ceitain standard deviation 
of the mean or probable range of erior. Thus, each 
computed valiieof Ehac! some uncertainty associated 
with it, The magnitude of the uncertainty depended 
both on the magnitude of the errois in each of the 
measured variables and on the nature of the equal ion 
itself. 

Because of the problems in evaluating D and Z 
accurately, the probable maximum erroi for each was 
assumed to be ± 25 percent of the absolute value for 
use in erroi analysis. The magnitude of these possible 
errors appeared reasonable in light of the preceding 
discussion. 

A probable uncertainty or ei ror in each calculated 
hourly value of E was evaluated by using procedures 
outlined in appendix A. All estimates of possible 
error were based on probable measurement errors for 
each variable (table 2). The estimated range or 
probable uncertainty in a particular E estimate 
assumed that all sensors and recording systems were 
operating within their accepted limits and that no 
error was introduced in the data-reduction 
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Table 1 ~A summary of conditions and measurements on days of experimental observations at 
lysimeter YW2C and watershed 109, Coshocton, Ohio 


Date 1 

1 

Time 

(EST) 

Soil 

condi' 

tions 

Accumu- 
Idted E 
(iysi- 
meter) 

Accumu- 

lated 

precipi- 

tion 

Cover 

and 

crop 

height (h) 

Av u 
(cm/s) 
direction 

General 
weather 
con tions 

Measure- j 
mcnts 

Height 
of Ae 
and aT 

i 

Length 

of 

fetch 




Cm 

Cm 

Cm 




Cm 

M 

5/2/69 

04'22 

Dry and 
loose 

on 

surface 

0 

0 

Bare 

soil, 

250 (SW) 

Clear and 

warm. 

Lysimeter 
Y102C,Rn. 
G, Ac, 

AT, All. 

to and 100 
D^O 

Zq = 0.1 

150 

5/22/69 

04-22 

Surface 

moist 

5.1 

6.1 

corn 

General 
bare soli 

corn 

h = 5. 

150 (N) 

Overcast Ly.sinictcr 
and cold Y102C,R„ 

0.05" rain G,Ae, 

@ 1200 EST. AT. Au. 

EST. 

10 and 100 

D = 0 
/.„ = 0 

15 

6/26/69 

04-22 

Wet soil 
saturated 
to field 
capacity. 

15.2 

20.9 

Corn 

almost 

complete 

cover 

h 75. 

200 (S\V) 

Clear and 

warm . 

Ly.slmelcr 
Y102C, R„, 
G, Ac, 

AT, All. 

150 and 250 

n = o 
*0 = 2.0 

150 

7/24/69 

04-22 

Wet 

26.7 

39.4 

Corn 
h = 240. 

100 (SW) 

Clear with 
haze, rain 
@1500 

EST 

Lysimeter 
YI02C,R„, 
G, Ac, 

AT, All. 

250 and 350 

D = 150 
z= 15.0 

150 

m7l69 

05-23 

Surface 
dry,. sub- 
surface 
moist. 

42.0 

52.0 

Corn 
h « 240. 

100 (N-S) 

Clear, 

light, 

variable 

wind. 

Lysimeter 
Y102C,Rn, 
G, AE, 

AT, All 
Wntcrslied 
109 AC, at. 

250 and 350 

D = 150 

Zq == 20.0 

15 to 150 

50 to 75 

10/14/69 

05-20 

Dry 


57.5 

Corn 

ready 

to 

harvest 
h= 130. 

150 (NW) 

Overcast 
and cool. 

Lysimeter 
Y102C, R„, 
G, Ac, 

AT, Til 
Watershed 
109, Rn 

Ac. at, 

All. 

150 and 250 

D = 50 

7.0 = 8,0 

15 to 

100 

10/16/69 

07-19 

Dry, 

surface 

moist. 


57.7 

Bare 

soil 

with 

litter. 

500 (SW) 

Partly 
cloudy 
and cool. 

Lysimeter 
Y102C,Rn, 
G, Ac, 

AT, Au 
watershed 

109, Rn, 

Ac, AT, 

Au. 

25 and 100 

D = o 
*0 = 0.1 

150 
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Table 2 — Recording anJ sefKsoi accur/icy of micronieteorohgfcai vafiables for use in probable-error analysis 


Variable 

Accuracy ot 

Probable 

incut 

Recorder 

Sensor 

error 

Rn 

±0 01 langley/iiiin 

±0 02 langley /mm 

±0 02 Janglcy/miti 

G 

± .01 langley /min 

± 02 langley /mm 

± .02 langley /min 


± 35° C 

± 10" C 

± .40" C 

AT 

± .02" C 

± .02° C 

± 03" C 

AT(wet bulb) 

± ,02" C 

± 02° C 

± .03" C 




± 60 mbar 

Ae 


- - - 

- .10 mbar 


— 

— 

^ .80 mbar 

1] 

±4 5 cni/s 

— 

'4.5 ern/s 

All 


— 

■*^6.4 cin/s 

Zq 

Assume 25% of 

-- - 

— 

D 

Assume 25% of 0 

... 



EXPERIMENTAL RESULTS 


Bowen Ratio and Aerodynamic Methods 

A comparison of hourly estimates E using the 
Bowen ratio and aerodynamic equations with the 
lysimeter value for May 2 is shown in figures 3 and 4, 
and for June 26 in figures 5 and 6» respectively. Both 
days were characterized by clear skies. Wind was 
from the south-southwest for most of the day» with 
speeds ranging from 100 to 300 centimeters per 
second. This direction provided the most ideal fetch 
conditions at the lysimeter test area. 

On May 2> measurements were made ovei bare, 
dry soil. As a result of the solar healing of the bare 
soil, measurements indicated strong temperature and 
humidity gradients, and the soil surface temperature 
reached 40° C during midday. Bowen ratios (B) 
ranged from L2 to 1.9 during daylight hours, 
indicating that most of the radiant energy was being 
transferred into sensible heat. 

The large negative Richardson numbers that 
prevailed during the measurements on May 2 
resulted in a large diabatic correction to the 
aerodynamic equation (8). Evaporation computed 
from the simple uncorrected logarithmic equation (4) 
averaged 50 percent less than those obtained with 
equation 8. This required the largest diabatic 
correction on any of the 7 test days. The closeness of 
the Bowen ratio estimate and the aerodynamic 
estimates using "‘KEYPS” correction was 
encouraging, There appeared to be a tendency for the 


aerodynamic equation with the Davis adjustment to 
overcompute E on this day. 

The range of probable erroi on May 2 was very 
small for both aerodynamic and Bowen ratio 
estimates, ranging from ± 0,01 to ± 0,02 langley per 
minute. The probable error for the Bowen ratio 
estimate was shown to increase considerably at 1800 
hours to a range of i 0,07 langley per minute. 

The run of June 26 was made over corn 75 
centimeters high following a week of heavy rains that 
had brought the soil to field capacity. This day was a 
marked contrast to conditions on May 2, in that the 
temperature profiles were very close to adiabatic 
during most daylight hours. Practically all the net 
radiation went into latent heat flux, with Bowen ratio 
values (B) ranging from -0.13 to 0.20 during day- 
light hours. Because of the low Richardson numbers, 
(-0,01 to +0.01), only a slight diabatic correction 
was necessary for the aerodynamic equal ion, ranging 
from 0.90 to I.O, maximum. 

Evaporation flux rates computed by aerodynamic 
and Bowen ratio methods on June 26 at the lysimeter 
site were compared with the measured lysimeter 
values (figs, 5 and 6), The Bowen ratio method 
agreed closely with the lysimeter data, while the aero- 
dynamic method provided slightly lower values 
during the hours 1 100 to 1600. 

Except for a couple of hours, the Bowen ratio 
estimate was generally within the probable error 
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LE (calories / cm^-min) 



Figure 3 — Hourly evaporalive flux from bare soil obtained on May 2, 1969, Irom weighinglysimetci Y [02C at Coshocton. Ohio 
and calculated with Bowen ratio method showing range ot probable error in estimating LH 


LYSIMETER SITE 



F^igurc4 —Hourly evaporative flux from bare sod obtained on May 2, 1 969. from weighing lysimeterY I02Cat Coshocton. Ohio, 
and calculated with adjusted aerodynamic methods sliowing range of probable error in estimating LE 
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LE (calories/cm —min) 



Figures —Hourly evaporative Dux from corn (75 cm) obtained on June 26, 1969, from weighing lysimetcr Y102C nt Coshocton, 
Ohio, and calculated with Bowen ratio method showing range of probable erjor in estimating LE. 




LE (caiories/cm^-min) 


LYSIMETER SITE 



Figure () — Hourly evnp oral i VC flux from corn (75cm)oblainccl on June 26, 1969, (lom weighing lysimctcr YI02C alCoshocion, OIiuk 
and calculated with adjuMcd netodynanuc methods showing lange ol probable error in cshmating LF 


range of 0.03 langley per minute from the lysimeter 
value. The range of probable error for the aero- 
dynamic estimate was very large from 1 100 to 1600 
hours, ranging from ^0.10 to ^0.15 langley per 
minute on an hourly basis. Reliability in the aero- 
dynamic estimate of E was gi eater during the early 
morning and late evening hours, 

D was assumed to be accurate within i 25 percent 
of the '*true” value for use in the error analyses. On 
June 26th, with h 75 centimeters and D = 50 
centimeters, this probable error in D would be* 12 
centimeters. If it were possible to estimate D without 
error, the range of probable uncertainty in the 
aerodynamic estimate of E would be reduced to 0.06 
langley per minute from 1 100 to 1600 hours. 

Thus, the large uncertainty in the aerodynamic 
estimate can be traced partly to the uncertainty in the 
estimate of D. Confidence could be improved by 


increasing the sampling heights so that D is small, 
relative to Z2 and Z i . However, if this is done, the 
profile measurements may not be representative of 
the processes occurring on the given surface. 

The comparisons described above are for the 2 test 
days considered most ideal. Conditions for the other 
days were often far from ideal. A summary of total 
daylight vapor flux (E) by the various methods for a 11 
periods is given in table 3 and is listed inappendix D 
alongwith the lysimeter values for the hourly periods. 

The Bowen ratio estimate was within 1 0 percent of 
the total lysimeter E measured on five of the seven 
cases and within 20 percent on the last 2 test days, The 
difference for the total E on these days is well within 
the accuracy of the meteorological and lysimeter 
measurements. 

The results obtained with the aerodynamic 
estimates were more variable. Only in a few cases did 
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Table 3,-A compamon of observed lysimeter E , with E computed hourly by Bowen ratio and 


aerodynamic methods with daylight total added 


Date 

Time 

(EST) 

E 

Lysimeter 

Y102C 

Bowen ratio 

Aerodynamic equation 

Fetch 

height 

ratio’ 

Equation 3 

Probable 

error 

Davis adj . 
equation 

8 and S 

Keyps adj 
equation 

8 and 7 

Probable 

error 

Keyps 



Mm 

Mm 

Mm 

Mm 

Mm 

Mm 





Lysimeter site 






5/2/69 

7-18 

1.42 

’1.40 

±0.08 

2.42 

*1.57 

--- 

150 

5/22/69 

12-20 

1.23 

M.33 

± .11 

*1 27 

’1 29 

--- 

15 

6/26/69 

5-22 

7.06 

*6.74 

± 12 

*6 35 

’5.90 

±0.37 

50 

7/24/69 

6-16 

3 81 

’3 46 

± 12 

*3.31 

2 30 

± .56 

45 

8/27/69 

6-23 

4.60 

’4 38 

± .16 

2.78 

1.96 

± .17 

5-45 

10/14/69 

8-18 

92 

.60 

± 08 

49 

.38 

± 04 

6 

10/16/69 

8-17 

1.10 

’ 90 

± .16 

*1.22 

*1.20 

± .01 

150 




Watershed site 





7/24/69 

6-16 

3.81 

’3 20 

± .11 

5.76 

*3.78 

±1.08 

20 

8/27/69 

6-23 

4 60 

’4 10 

± .21 

3.05 

2.10 

± 37 

15-20 

10/14/69 

8-18 

.92 

’ 94 

± 17 

M 30 

M.IO 

± .08 

40 

10/16/69 

8-17 

1.10 

.46 

± .5 

08 

.07 

± > • • 

75 


^ Ratio of minimiim upwind distance from sampling site of surface being measured to height of highest sampling point in profile 
measurements used m computations. 

’ Computed E within ±10 percent of lysimeter E. 

^ Computed E within ±20 percent of lysimeter E 


the computed flux rates closely agree with the 
lysimeter or Bowen ratio estimates of E after crop 
height had exceeded 75 centimeters* Between May 2 
and July 24, the probable range of uncertainty 
increased from * 0*07 to ±0.56 langley per minute, as 
D increased because of crop growth from 0 to 150 
centimeters. The total flux estimate on July 24 would 
have had an uncertainty of ±0.24 langley per minute 
if it could be assumed that “D” had been measured 
without error. Thus, a major source of error in the 
aerodynamic method appears to be the uncertainty in 
the measurement or approximation of the displace- 
ment height D. 

A large variation occurred in the aerodynamicand 
Bowen ratio estimates at the lysimeter site on 
August 27 and again on October 14. This failure can 
be related to poor fetch conditions; light, variable 
winds arose on both days. On October 1 4, the water- 
shed site was better suited in regard to wind direction 
and fetch. The results at the watershed site showed 
that all three estimates agreed on this day. In 
contrast, on October 16 and under heavy winds, the 
watershed site failed to provide any meaningful 
estimates, whereas the lysimeter site with a longer 
fetch closely agreed on all three estimates. Generally, 
the Bowen ratio estimates of E provided reasonable 


accuracy when the fetch / height ratio was in the range 
of 20 to 1 through 50 to I, except for October 16 
when, for bare soil, a fetch of 75 to 1 was inadequate. 
A fetch/height ratio of greater than 50 to 1 appears 
necessary for successful application of aerodynamic 
techniques. 

Combination Formula 

Four days— May 22, June 26, July 24, and 
August 27 — were considered to have met the require- 
ments of nonlimiting water supply necessary to 
satisfy the definition of potential evapotranspiration, 
Latent heat flux was computed hourly for each day 
from equation 9 using Zq as detei mined from wind- 
profile analysis. The results for two of these days are 
plotted in figures 7 and 8. On June 26, the estimated 
value of Ep was close to E, but about 10 percent 
greater than the amount as measured by the 
lysimeter. The range of probable error in the Ep was 
generally ±0.04 to ±0.08 langley per minute, 
assuming an uncertainty in Zo of 0.5 centimeter, 
with Zq equal to 2,0 centimeters. Thus considering 
the range of uncertainty, hourly Ep estimates agreed 
with the lysimeter except, perhaps, from 1400 
to 1700. 
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LYSIMETER SITE 
JUNE 26, 1969 



TIME END OF PERIOD (hrs. E.S.T) 

pigurc? —Hourly evaporative flux from corn (75 cm) obtained on func26, (rom weighing lyhinicler Y)02CatCc>shoclon, 
Ohio, and calcidaled potential El" by combinaiion cijuation showing range of probable error irt estinumnig 1 I 

LYSIMETER SITE 



TIME END OF PERIOD (hrs E.S.T) 


Figures — Hourly evapoi alive flux from corn (240 cm) obtained on August 27, 1969, from weigh in glysimc ter Y I t)2C a I Coshoc- 
ton, Ohio, and calculated potential E'l by combination equation showing range of probable error in esiimahng I I* 
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In contrast, the lesults foi August 27 revealed an 
overestimation of E by moie than 100 percent of the 
I ysi meter value. Ep also exceeded net radiation by 10 
to 15 percent fj om 1 lOO to 1600. The range of 
probable erroi (±0.15) was very large compared with 
the June 26 range in probable eiror (±0.06). 

A summary of totalized daylight hourly values for 
each of the 4 days is given in table 4, column 2. 
Equation 9, using measured values, failed to 
predict E accurately once the corn had reached full 
height at 240 centimeters. On May 22, agreement 
between computed Ep and lysimeter flux rales for 
bare soil was excellent despite poor fetch conditions. 

The value can be estimated as a function of crop 
height (26) by the equation. 

= 0.058 h‘-19 

The value Z (j could have been estimated on the basis 
of equation 1 1 rather than on the basis of a 
transposed value from the Ithica, N.Y,, experiment 
previously mentioned or wind-profile analysis. Then, 
the values of Z o would be 10 and 40 centimeters for 
crop heights of 75 and 240 centimeters, respectively. 
Used in equation 9 these Zo values would greatly 
increase the difference between computed E p and 
observed E. Thus, it was necessary to substitute the 
value from the Ithaca, N.Y., expeiiment because of 
the nonhomogeneous terrain at the Coshocton, Ohio, 
site. 

In reference to the totalized daylight values given in 
table 4, the probable error in the period total for 
individual days increased considerably as rose 


from 0. 1 to 20 centimeters, assuming an uncei tainty 
of ±25 percent in the Zq estimate. An absolute 
measure of Zo would reduce the range ofpiobabihty 
error to ± 0. 1 3 millimeter per day or less for all 4 days. 
This uncertainty in the E p estimate could account for 
some of the difference; howeveu it is doubtless not 
the only reason for the overcomputation of E on 
July 24 and August 27. 

The log law for neutral profiles fiom which 
equation 10 was derived is recognized to be in error 
for diabatic situations. Tanner {26) suggests that the 
wind function (10) includes a diabatic convection 
factor similar to that used in the aerodynamic 
equation coirection. However, if such a correction 
factor were applied to the 3 days in question, it would 
increase the error between estimated and observed 
lysimeter flux lates. Accounting for this disagree- 
ment would be easy by assuming that the surface was 
not a potential evaporation surface or that the 
measurements of the meteorological variables were 
not representative. The success of the Bowen ratio 
estimates invalidates this later argument. 

If the surface is not evaporating at the potential 
rate, then computing the surface diffusion resistance 
(i which is a means of defining the degree of 
unsaturation of the crop and soil system, should be 
possible. When the energy balance of the crop and the 
aerodynamic resistance term are known, rs can be 
calculated as a lesidual from the ratio of potential to 
actual evaporation (25). 
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Table 4. -Observed lysimeter E during periods of adequate moisture supply, compared with potential 
evaporation (Ep), computed hourly by combination equation with total for period 


Date 

Time 

(EST) j 

E 

Lysimeter 

Y102C 

(1) 

Equations 

9 and 10 
measured 

(2) 

Probable 

error 

Zq ±25% 

(3) 

Equations 

9 and 10 
assume 

Zq = l .0 cm 

(4) 

Penman 

version 

(5) 

Measured 

^0 

(6) 



Mm 

Mm 

Mm 

Mm 

Mm 

Cm 

s/22/69 

1 200 to 2000 

1.23 

'1.06 

±0.07 



'1 17 

0,1 

6/26/69 

500 to 2200 

7.06 

’7,85 

± .20 

‘7.14 

'6.78 

2 

7/24/69 

600 to 1600 

3.81 

4 55 

t ,28 

'3.51 

'3.60 

15 

8/27/69 

600 to 2300 

4.60 

8.14 

t .47 

'4.87 

'5 03 

20 


‘ CompiUed Ep within ± 10 percent of lysimeter E. 
^Computed Ep within ±20 percent of lysimeter E. 
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Potential evaporation (Ep ) can be calculated from 
equation 9, When the measured Z o values were used, 
the calculation of r., provided hourly values from 0 to 
0.8 second per centimeter for June 26, and from 0.5 
to 2.5 seconds pei centimeter for August 27 during 
daylight periods. These estimates for r*, appear 
reasonable, but there are no independent measures of 
r^ orsoilwatei potential foi comparison Because of 
the high degree of lecharge during this period 
(table 2), precipitation - measured evaporation), 
probably no soil moisture deficiency existed. Thus, 
water was not limiting, and potential conditions 
should be as valid on July 24 and August 27 as on 
June 26. 

Jensen (9) and Rosenberg (27) both report failure 
of the combination equation 9 when using the 
aerodynamic term (10) proposed by Van Bavel ( 31 ) 
Jensen indicates that under arid conditions it is 
necessary to calibrate r., by vaiying the Zq value 
until the Ep estimate of equation 9 had agreed with a 
lysimeter measuie. Sellers (22) has indicated that 
unreasonably high E estimates have been obtained by 
using simple logarithmic profile methods over tall 
ciops. Van BavePs verification of the modified com- 
bination equation (57) was limited to bare soil and 
alfalfa 25 to 30 centimeteis high. The authors too 


could find no reported applicationsof equation 9 ton 
tall crop such as corn, where an independent measure 
of E was available for verification. 

Potential evaporation was recalculated for the 3 
days with cropcovei byequation 9, usingan assumed 
Z„ of I.O centimeter in the wind-function term (10). 
Improvement in the period totals for each of the 3 
days can be observed in table 4, column 4. While the 
change in Zo from 2 centimeteis to I for June 26 
resulted in a slight improvement, the reduction in Zq 
from 20 centimeters to I on August 26 resulted in an 
Ep estimate that is within 5 percent of the lysimeter 
v^ue, as compared with a pievious ovei estimate of 
100 percent. These results agiee favoiably with those 
of Jensen (9), who reported that for a variety of crops 
in southern Idaho, in equation 10 should be 
approximately 0.75 centimeter where crops were 
higher than 50 centimeters. 

For compaiison, the combination equation 9 
using Penman’s original aerodynamic term (77) was 
computed hourly and summed to give a total for each 
day. The results are shown in table 4 column 5. They 
indicate that the wind function as originally used by 
Penman produces essentially the same results as the 
log profile version (10) using a Zq value of 1 
centimeter. 


DISCUSSION 


Most of the meaningful verifications of the Bowen 
ratio energy budget technique have been reported for 
low-growing inigated crops with generally flat 
topography and unifoim fetch conditions. The 
success of the hourly values of E computed by 
equation 3 for the lysimeter at Coshocton, Ohio, 
tends to verify the application of the approach to 
estimating evapotranspiiation under natural field 
conditions where fetch and surface homogeneity are 
less than ideal. The method is capable of predicting 
actual evaporation from a bare, dry soil with limiting 
moisture, as well as from a tall ciop such as corn 
(height “ 240 centimeters) under a wide range of 
growth and moisture conditions. 

Overcast skys (low direct solar ladiation) or 
variable weather conditions do not appear to limit the 
application of this technique except, perhaps, during 
periods of rain. The Bowen ratio method fails during 
early morning and late evening as a result of 
increasing relative errors in measuring temperature 
and humidity profiles during nondaylight hours. 
However, since E rates are quite small and are often 
zero at such times, this failure is relatively 
unimportant. 


Use of aveiage houily values of AT and Ae, as 
measured by the Lourence and Pruitt psychromctcr 
system ( 14 ), appears sufficient to give accurate 
estimates of B for the Bowen ratio technique. The 
success of the Bowen ratio method for the estimation 
of hourly E flux rales indicates that the micro- 
meteorological measurements of radiation, temper- 
ature, and humidity are representative of conditions 
prevailing at the test site. 

The success of the Bowen ratio methods under 
poor fetch and sometimes variable conditions 
indicates that fetch requirements are less strigent for 
this approach as compared with the aerodynamic 
methods. To give accurate results, the required height 
of the sensor-to-fetch-distance ratio appears to be no 
greater than 1 to 50 and possibly smaller. On two 
occasions, one over corn and the other over bare soil, 
a height-fetch ratio of about 1 to 10 was adequate. 
However, on October 16, with wind speeds greater 
than 500 centimeters per second, a ratio of 1 to 75 was 
inadequate at one site, whereas the other sampling 
site with a ratio of 1 to 150 provided an excellent fit 
with evaporation measured by the lysimeter. 
Perhaps, during this last sampling period a greater 
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sampling frequency and a shorter computational 
period would have impioved the estimate of E by 
equation 3 

On the basis of results with aerodynamic estimates 
of E fetch distance appeared to be fiom 75 to 100 
meters for light and moderate wind conditions, up to 
J50 meleis or greater for winds in excess of 500 
centimeters per second when the upper piofile 
measurement height was 100 centimeters above the 
crop surface A diabatic correction function must be 
applied to the accepted Thornthwaite and Holzman 
aerodynamic expression (29), equation 4, duiing 
no nadiabatic conditions in order to agree closely with 
lysimeter and Bowen ratio E rates. The Keyps and 
Davis empirical correction functions appear to be 
equally satisfactory; however, occasionally a 
significant difference between the two approaches 
occurs. No recommendations can be made on the 
basis of the tests described here as to the advantages 
of one correction function over the other. Generally, 
the aerodynamic method gives better estimates of E 
than does the Bowen ratio method only during eaily 
morning and late evening hours. 

The error analysis tends to indicate that much 
uncertainly in the aei odynamic methods results from 
error in measuring or in estimating D. This eiioi is 
shown to increase in range as both D and crop height 
(h) increases. Thus, unless D is carefully determined 
by accurate wind-profile analysis, the aerodynamic 
methods may have a large error when applied ovei 
tall vegetation. Thus, aerodynamic methods appear 
to require better conditions 

Estimates of evapotranspiration with the improved 
versions of the combination equation, using 


value.s based on wind piofilc oi crop height, are 
unsatisfactoiy once ciop height exceeds 75 
centimeters A 25 peiceni unceitainty in icsults 
in a probable ei 101 of ±20 percent in E (August 27), 
whenZ(> = 20centimetcjsascompaied vvitha 5 per- 
cent error in E if Z^^ weie assumed to be measured 
without enor. Since the unceitainty in the 
determination of Z^ increases with crop height, 
resulting m a high degree ol unceitainty m potential 
evapotranspiration estimates over tall vegetation, the 
use of the log-profile wind function should be limited 
to bare soil and to crops below 75 centimcteis on Hat 
terrain, 

Empirical methods can be used to fit a us liable 
wind function in the combination equation. IfZ^^ in 
equation 1 0 is used as a calibration tei m as suggested 
by Jensen (9), then Z^, in this application should be 
1.0 centimeter foi corn as high as 240 centimeters. 
This proposition assumes that no icsistance to vapor 
transpoit within the corn plant occur icd and that the 
potential evapotranspiration concept is valid as 
defined. The empiiical coefficients proposed by 
Penman foi the wind function of the combination 
equation appear to be generally applicable to a tall 
crop such as coin in a humid region. Penman’s 
veision (J7) gives essentially identical icsults to Van 
Bavel’s veision (31), using a Z^ of 1.0 centimctei. 

Two conclusions aie possible from the limited 
results with the combination equations: ( I ) There is a 
resistance to water vapor transport within the coin 
plant despite high soil moisture conditions or (2) the 
log-piofile wind function is unsatisfactoiy for 
application in the combination equation whcie crop 
height exceeds 75 centimeters on heterogeneous 
terrain. 
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LIST OF SYMBOLS 


Symbol 

Explanation 

Units 

B 

Bowen ratio, B = H/E 

1 

D 

Displacement parameter, D = d + 

cin 

E 

Evaporation energy equivalent 

2 1 
cal enr^s"^ 

"p 

Potential evaporation 

1 “2-1 
cal cm s 

G 

Heat flux into ground 

1 “2 -1 
cal cm s 

H 

Heat flux into an 

2 1 

cal cm s 

Kh 

Eddy diffusivity for heat 

2 -1 
cm s 

Km 

Eddy diffusivity for momentum 

2 -1 
enr s 

Kv 

Eddy diffusivity for vapor 

2 “1 
ciir s 

L 

Latent heat of vaporization 

585 cal/g 

P 

Atmospheric pressure 

mbar 

Kn 

Net radiation flux 

cal cm s ^ 

Tz 

Temperature at height (z) 


Z 

Height above ground surface 

cm 


Specific heat at constant pressure 

cal giTi“^ °C' 

cl 

Zeio plane displacement 

cm 


Actual vapor pressure of air 

mbar 

«s 

Vapor pressure of air at satuiation 

mbar 


Vapor pressure of air at height (z) 

mbar 

g 

Acceleiation of gravity 

cm 

h 

Vegetation height 

cm 

k 

Karman’s constant 

1 

In 

Natural logarithm 



Aerodynamic diffusion resistance 

s cm'* 
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LIST OF SYMBOLS-Continued 


Symbol 

Explanation 

Units 


Surface diffusion resistance 

s enr^ 

s 

Slope of saturation vapor pressure 
curve s = de/dT 

mbar 

“z 

Horizontal wind velocity at height (z) 

cm 


Roughness lengtii 

cm 

5 

Wet-bulb depression 

"C, °K 

e 

Ratio molecular weight water vapor to air 

1 

7 

Psychrometcr constant, == CpP/Le 

mbai C“' 

P 

Density of moist air 

-3 

g cm ^ 


Stability profile influence function 

I 

A 

Increment 


Subscripts 




Property at height (z) referenced above ground level. 
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APPENDIX A-PROPAGATION OF ERRORS 


To determine the nature of error propagation ( 4 ), 
as error analysis is sometimes called in indirect 
measurement, consider the quantity Y which is a 
function of independent variables , X 2 , X n 
(Y = f[Xj,X 2 . . ^n]) The change in Y that is 

due solely to a change in X 1 may be denoted by 
3Y/dX| , where dY/5Xj is the partial derivative of Y 
with respect to Xj and where dY/^Xj represents 
the change in Y per unit change in X| , while dxj 
lepresents the actual small change in X| The total 
change in Y that is due to changes in all the variables 
is denoted by the symbol dY,n , where: 


= 

m 


9 Y 


0 Y , 


aY 

— dx. 

+ 

— cix. 

. . + 

— dx 

9 Y ' i 




n 


(lA) 

Thus, if dx| , dx 2 , . . . dxn represent the 
uncertainties in the directly measured quantities X| , 
Xj,. . x„, respectively, and if all the uncertainties 

aie considered additive, then dY„, represents the 
maximum uncertainty in Y. This uncertainty in Y 
may be expressed in terms of any desired measure of 
dispersion, as long as the same measureofdispeision 
is used for the uncertainties in the diicctly measured 
quantities also, Howevei, it is extremely unlikely that 
the actual error in a calculated quantity will eciiial the 
maximum uncertainty expressed by equation I A 

The uncertainty or possible erroi in a calculated 
quantity is repiesented better by a probable value 
than by a maximum value. The probable uncci tainty 
(dYp ) or error { 4 ) in a calculated quantity (Y) is 
determined by the square root of the sum of the 
squares of the individual terms in relationship to the 
maximum uncertainty. The resulting equation is; 


dYp = 


av 2 2 3Y 2 ,, 

dX|^ + (- — )2 dXn^-t,, + 
dX| 9x2 


9 Y ~ , 

n 


1/2 


(2A) 


Thus, the probable uncertainty or error in Y 
depends not only on the uncertainty of each directly 


measured quantity but also on the partial derivatives 
of the equation that repiesents Y. A vai iable would 
not contribute towaid the probable error in Y, 
regardless of the value of dY/5xi,ifthc uncertainty 
(dx ] ) for that variable were zero. The nature of the 
equation docs not contribute directly to the erroi but 
rather modifies the enots already present in the 
measured values. 

The following example details the proceduie used 
to evaluate the probable and the maximum error in 
the Bowen ratio estimate (equation 3) of E 

By taking the partial derivatives of equation 3 
successively with respect to G,AT, andAe, and 
by substituting the i esiilting expressions in I A and 
2A, one can obtain the piobable erroi: 



l+B 


1^ dR.,“) + (I 


1 +B 


1 “ cIG") + 


(- Rq + G> B ^dAT^ 

( — (—) ) + 

L(l + U)2 AT J 


( 


(Rn+G) 

(1+B)2 



dAe^ ) 


{3A) 


and foi the maximum uncertainty 


dE„ = 


-L) cIR,, 
1+B " 


,(_L) dG 
' 1 +B 


-(Rn+G) B 


(Rn+G) B 

tlAT 

+ 

-) dAe 

(l+B)'^ AT 


(l + B)-^ Ae 


(4A) 


where , dG, dA'i , and d represent the piob”'^'" 
error in measurement of the variables R„ , G, A 
A,e, respectively. 



On May 2 , between 1300 and 1400, the following 
values were observed: 


R ^ 0 620 cal/cnrmin. 
n 9 

G = -0 055 cal/cm - mm 


AT,„.,00. ■ 


Ae 


1 0-1 00. 


1.507 mbar 


E = 0 221 cal/cm . - mm 


The value of the probable measurement error for 
each variable was derived from table 2. By 
substituting this error in equations 3A and 4A,one 
derives the following error estimates. 


dEp - i 0.012 cal cm - mm 
dE^ = ± 0.018 cdlcm^ - mm. 


Thus. E 3 1400 = 0.221 ± 0.012 cal/cm^ -min. 
Equations similar to 3 A and4A were developed from 
equati ons 8 and 9 by the procedure outlined above for 
error analysis of these estimates. 


The maximum unceitainty of E is: 


dHn, = ± 



— 

A Au 





'^2-D 


2 d Ae 


i 

1 

D 
1 

) 



1- \ 





A Ae 



dAu 

L('tvX 3 j)J 



2 A Ae Au (^2 - Z | 


dD 


L ■% 


3 (Zj -D) (Zj -D) 


(7A) 


Stability Correction 


Aerodynamic Equation 
E = Ae Ati 

( Jn (Zj-D/Zj-D))^ (5A) 

Letting A = eXp and assumingthatZ is measured 
without error, the probable error in estimate of Y E 
is; 




“ 

r A Au -| 

2 


\Jh -D\2 

2 


L \"! - 

(lAe 





AAe 


In 


Z2-D 


Zj -D 



clAu 




1 - 


bgAT AZ 




(8A) 


The probable error (d foi the estimate of i// is; 


ci-//p = ± 


bg AZ 


4 .,.5 




2 iIAT^ + 


bgAT 


(Au)2 


dAZ'^ + 


AT AZ 




"^abs 



r 2 A Ac Au (Z2-Z() 

2 2' 

dD 


r "1 


\n (^2 ■®\(Zt -D)(Z, -D) I 


(6A) 1 

-bg AT AZ 


- J 



2^^ (Tabs) (A«)^ 


d Au-* 


(9A) 
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The maximum uncertainty of \p is: 


fj 0 i hg 






4 1 //“ 


dAZ 


AZ 




cIT 


AT AZ 


cIT, 




abs 


-2AT AZ 

JabsA')^ 


d All 


(lOA) 


r 2pCpK2(e^-e^,)„ 2 


LS 


Z-D\13 


In 





The maximum uncertainty of E p is: 


‘*(Ep),n 


, S 
S+7 


dR,. 


ilG 


(I2A) 


Combination Equation for Potential ET 


= 


S+7 


pC,,k u(c„ -c„) 
(R„ + G) + _JL ^ 


2 

(IIA) 


Assuming that S can be calculated without error, the 
probable error in estimate of (E p) is: 


" S+7 


d R„2 (It; 2 





du 


-h 





I" 0CpK2(o^.„p 


LS 


In/ Z-dV 2 


2di|2 ^ 


1 

to 



-2 " 

2 d(C3-e_^)2 

+ 

LS 


'■0 . 





2pCpK2(c^ -Cj,)!! 


LS In 


Z-D 3(Z - D) 


2 

do2 


+ 








(13A) 
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APPENDIX B-INSTRUMENTATION 
AND DATA REDUCTION 


Lysimeter 

The weighing lysimeter that was used to measure E 
directly has already been described m detail (7, S) 
The wide concrete perimeter between the lysimeter 
and the surrounding field, as originally constructed, 
has since been removed, and crop planting has been 
brought to the very edge of the lysimeter. The 
automatic weighing mechanism records weight 
changes in the 59,000 kilogram mass to a 2.26- 
kilogram accuracy, which is equivalent to 
approximately 0.25 millimeter of water ovei the 
lysimeter surface area, Average hourly E was 
obtained by averaging the six consecutive 1 0-minute 
weight s begi nni ng at 60 mi antes bef ore each li ou r a nd 
ending at the hour. This averaging process in weight 
records removed wind-caused irregularities. 

Wind Profile Measurement 
and Data Reduction 

The wind profile measurements at the lysimeter site 
were made with four matched rotating cup 
anemometers that had been manufactured by Wong 
Laboratories. These were mounted 50, 100, 150, and 
200 centimeters above the crop surface. The sensors 
produced electrical contact closure for every 0.1 mile 
of wind passage, sending an impulse to a Sodeco 
printing counter. The total number of impulses for 
each hour was printed on the hour, providing the 
average wind speed in miles per hour. 

Wind-profile measurements at the central 
watershed site were made during three of the seven 
runs. During the July 24 run, a system of six 
Thornthwaite rotating-cup anemometers was used. 
Impulses were counted on mechanical counters. 
These anemometers produce an impulse for 
approximately every 1,5 meters of wind movement. 
During the runs of October 14 and 16, a Beckman 
and Whitley system of three rotating-cup 
anemometers was used. Impulses were counted on 
mechanical counters that were photographed on 
microfilm every 15 minutes. These anemometers 
produced an impulse for approximately every 1,5 
meters of wind movement. The total number of 
counts for each hour was recorded and converted to 
average wind speed in centimeters per second. 

The methods used to evaluate D and were 
subject to shortcomings and required sensitive 
instrumentation, suitable fetch, and neutral 


atmospheric conditions. Temperature gradients and 
elastic properties of vegetation combined to make D 
and z difficult to cliaracteiize for vegetated surfaces. 
Since temperature gradients are relatively small over 
vegetation, plant geometric properties are the 
greatest influence on D and Therefoie, the 
parameters can be expressed as a function of crop 
height as shown in figure 9 (/2), 

d (cm) 



higiire 9.— Roughness length (zq) and 7ero plane displace- 
ment (D) as a function of height of ihe corn ciop (h) 
during the 1961 glowing season at Ellis Hollow, N Y (/2), 
Values aie for near-isotlicrmal conditions and wind speeds ol 
l(X) to 400 cm/s at the l-mctci level above the ciop siu lace 

Wind-profile parameters were determined by 
graphical data where In (Z-D) is visually plotted 
against wind speed and forced through intercept 
Zq values. Since the data points were often too 
scattered to provide a precise determination of Zo 
and D, these values were first estimated on the basis 
of wind-profile parameters measured over corn of 
equal height (h) at Ithaca, N.Y, (figure 9). This 
transformation is possible because of the sim ilarity in 
'‘Leaf Area Index” (LAI) between the corn at Ellis 
Hollow, N.Y,, and the corn used in this study 
(table 5). Except for I day (October 14), the 
transposed D and Zq values provided a good fit to 
match the hourly wind-profile data. On October 14, 
the D and Zothat gave the best straight line fit 
between In (Z-D) and wind speed were chosen. 
Values of U 2 andU| were taken from the straight line 
plot for each hour of heights (Z 2 - D) and (Zj - D) 
for use in appropriate equations, 


24 




Table 5 -A comparison of leaf area hides (LAI) at 
Ellis Hollow, N, Y,, and Coshocton, Ohio, locations 
as a function of relative height (Zjh) 


z 

Z 1 

li 

S LAi 

Clhs Hollow' 

LAI 

Coshocton* 

240 

1 00 

0 

0 

227 

.95 

.08 

.08 

189 

.78 

.40 

.40 

135 

56 

1 90 

1.90 

81 

34 

3 20 

3.20 

37 

15 

4 00 

3 60 

0 

0 

C) 

C) 


' Accumiildted LAI from top of crop downward 
^ Base level of Ellis Hollow, N Y , corn crop - 4 5 LAI. 

^Base level of Coshocton, Ohio, corn crop; LAI not known 


Wet" and Dry-Bulb Air Temperature 
Instrumentation 

An aspirated psychrometer system modeled after a 
system developed and tested at Davis, Calif., {14) was 
used to measure temperature and humidity profiles. 
The wet-bulb air temperature measurement and (he 
water supply were redesigned from the Davis system. 
Changes weie also made in the mast and clamping 
system, and minor modifications were made in the 
sensor-holding assembly Two separate profile 
systems were used: (1) A 10-level profile system in 
midwatershed site and (2) a 5-tevel profile system at 
the lysimeter site. In the five level piofile system, all 
measurements were taken above the crop canopy. 
The psychrometer systems provided measiiiements 
of temperature and humidity accurate to ±0.020*^ C 
and ± 0, 10 millibar, respectively. 

Figure 1 0 presents a cutaway di agram o f a di y- a nd 
wet-bulb assembly. The dry-bulb thermopile wires 
are encased within a thin-walled, stainless-steel tube, 
five-sixteenths of an inch in diameter, These wires 
emerge from the upwind and extend 4 centimeters 
beyond the end of the tube. Because the cool wet-bulb 
element lies parallel to this tube for some distance, 
some cooling effects can lesult.’ Having a tube of low 
thermal conductivity minimizes the rate at which heat 
can be conducted along this tube from the region 
where the dry-bulb thermopiles are located, A rubber 
sleeve shrouds the dry-bulb thermopile elements that 
extend from the stainless-steel tube. This rubber 
sleeve also fits snugly over the end of the stainless- 
steel encasement tube. The wiies are encapsulated 
within this sleeve with silicone rubber. 

The dry-bulb sensor consists of a five-junction 
thermopile set referenced to the dry-bulb element of 



higurc 10— Cutaway view of wet-bulb iincl dry-bulb Jbscmbly 
used in psychrometer system 


Other sensois. An additional independent feve- 
Junction thetmopile set extends back through the 
mounting plug and into the wet-bulb element con- 
tained on the same sensor plug assembly. The dry- 
bulb temperature differential system shown in 
figure 1 1 permitted a more precise measurement of 
the dry-bulb temperature profile and use of the 
recording system already on hand. At the lowei level 



Figure I'l. — Dry-bulb thermopile pro! 



sensor, dry-bulb temperature was referenced to an ice 
bath located at the base of the mast, or thediy-buib 
temperature was referenced to ice with a small 
thermocouple that had been encased along with the 
dry-bulb thermopile 

The wet-bulb element, which is shown in figure 1 2, 
IS not a conventional cotton-wick type, but a porous 
ceramic cup 14.6 centimeters long. 

Water is supplied to the porous cup through a clear 
plastic tube that leads from a supply bottle located on 
the main sensoi;-mounting arm. This tube is 
connected to a 0.79-cenlimeter (five-sixteenth of an 


inch) diameter stainless-steel tube that passes 
through the sensoi assembly mounting plug. The 
porous wet-bulb element is attached to this stainless- 
steel tube with a threaded-brass connectoi . 

One advantage of the ceramic-cup, wet-bulb 
element is a closed water system. This system 
eliminates part of the problem associated with the wet 
cotton-wick system, wheie a fiee water surface was 
required to supply moistuie foi evapoiation on the 
wick The ceramic wick does not flood when the 
elevated waterhead is subject to evaporation, and 
when the wateihead is slightly negative, watei still 




END VIEW 

Figure 12 —Details of censor unit and aspirated shield used on psych rome ter system. 
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continues to move into the ceramic element because ' 
the element acts as a tensiometer with a continuous 
water column. However, one must bleed all the 
entrapped air in the wet-bulb water system before 
starting to opeiate. Under most field conditions, the 
water bottle reservoir does not require refilling more 
often than every 24 hours 

The wet-bulb thermopile wires run inside the 
porous ceramic cup to a point near the forward end 
where the sensingjunctions are located. These wires 
leave the water-filled media through a hole drilled at 
an angle from the back of the epoxy mounting plug 
and through the side of the 0.79-centimeter (five- 
sixteenth of an inch) diameter, stainless-steel, water- 
feed tube. This exit hole is later sealed by refilling 
with epoxy after the sensor wires are pioperly 
positioned. The five thermopile wires extend into the 
dry-bulb lube that has been described earlier. The 
two wet-bulb thermopile output lead wires originate 
from the wet-bulb end and pass from the sensing plug 
with the thermopiles and follow the dry-bulb 
thermopile wires that lead diiectly to the recording 
system. 

Figure 12 shows a scale drawing of a sensor unit 
and its radiation shielding tubes. A double shield 
arrangement that has been patterned after the Davis 
design is used. The outer shield consists of an epoxy 
fiberglass tube 0.12 centimeter (0,047 inch) thick, 
while the inner tube is a standard 1 -inch-type I PVC 
pipe. The outer surfaces of each shield are covered 
with a layer of adhesive mylar-coated aluminum foil. 

Air is drawn past the dry- and wet-bulb sensors at a 
velocity of 500 centimeters per second. A small 1 15- 
VAC fan capable of moving 1 5 cubic feet of air per 
minute while requiring 7 watts of power does this. Air 
is also drawn between the inner and outer shields at a 
rate of 16 centimeters per second through six small 
bleed holes located on the back end of the inner shield 
tube. This aspiration insures removal of any heat 
loading caused by radiative or conductive sources. 
Figuie 13 shows complete dry- and wet-bulb unit 
along with all of its disassembled components. 

A test was performed to determine the time- 
response characteristics of the dry- and wet-bulb 
elements in order to achieve good phase relationship 
between the wet-bulb and the dry-bulb temperature 
responses when the units were subjected to abrupt 
changes in temperature. A time constant of 20 to 25 
seconds after an abrupt change of 7,5® C was 
necessary to achieve 62 percent of the final 
temperature. During construction, each thermopile 
junction was checked for proper insulation and for 
verification that its response at a given temperature 
was within ± 10 microvolts of one another. The final 





Figure 13.— Complete dry- and wet-bulb iinii along with Us 

disassembled components 

calibration constant for the 5 junction thermopile 
assemblies was: 

°C « 4 876® C/mv 

Two profile masts were outfitted with the 
temperature sensors described. Figure 14 shows the 
level masts with the various sensor spacing used, A 
Condorf steel-channel tower 1 0 meters long was out- 
fitted with sensor units, The sensor units were 
mounted with I -inch clamps, AUenwrench set screws 
permitted adjustment of the sensor units in a vertical 
direction. The dry-bulb reference thermopile sets 
extended about 2 meters from the mast’s base, which 
had been placed in an ice reference bath contained in 
a large, foam-plastic ice pail. 

The differential millivolt signals for temperature 
difference and wel-biilb depression were recorded on 
a multichannel millivolt recorder (L &N Speedomax 
Type G), A summary of these measurements is 
contained in table 6. The recorder was connected by 
shielded-copper leads to the differential thermopile 
system. 

Instantaneous microvolt readings of temperature 
difference and wet-bulb depression were taken from 
the strip chart records at 5-minute intervals and 
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Figure 14 —Complete five- level, dry- and wet-bulb temperature 
mast. 


averaged on an hourly basis. These data were then 
converted from microvolts to degrees Celeuis, and 
average hourly dry-bulb and wet-bulb depression 
temperatures were computed. By knowing the dry- 
bulb temperature and wet-bulb depression at each 
sampling level, computing vapor pressure of the air 
for each hour period was then possible. 

The dry-bulb air temperature (°C) and vapor 
pressure (g/m^) were plotted against the In (Z-D), 
where D was determined from wind-profile analysis. 
A straight line was drawn visually to provide the best 
fit for the data points for a given hour. Values of € 2 , 
e I ,T 2 , and T^ were taken from the straight line plot 
for each hour at heights (Z 2 - D) and (z^ - D) for use 
in the required computations. 


Description of Other Temperature 
Measurements and Data Reduction 

Other temperature measurements made at 
lysimeter Yioic are summarized in table 6. Briefly, 
all absolute temperature data were recorded on a 16- 
channei Leeds and Northrup speedomax G 
thermocouple recorder with a range of -20® to 120^^ F 
and a sensitivity of 0.25® F. Dry-bulb and wet-bulb 


temperature measurements were made in each of the 
five aspiiated sensor assemblies described in the 
temperature and humidity piofile section. These 
measurements piovided an absolute temperature 
measurement foi the theimopile psychromelic 
system and seived as a check in comparison with 
absolute tempera tui e comput ed for each profile level . 

Instantaneous tempciature data were taken from 
the strip chart records at lO-minute intervals and 
tabulated for each 24-hour period. All readings were 
conveited fi cm degrees Fahienheit to degrees Celsius 
and average hourly values were computed fiom the 
10-minute values. 

Soil surface temperature was measuied by the 
Barnes IT3 infrared thermometer. The infraied 
thermometei instrument calibration was obtained by 
observing the millivolt output of the theimometer 
sensor while reading a black-body surface of known 
temperature. The complete curve for the working 
range of the instrument was detei mined by lepeating 
the procedure at various tempei atures in a water bath 
calibration source that duplicated black-body con- 
ditions. The output signal fiom the IT3 console was 
recorded on a Royson electronic integrator that had 
an accuracy of i 0.1® C. for the complete system The 
millivolt output from the Barnes infrared 
thermometer was averaged for each hour, and the 
equivalent temperature in ®C was determined from 
the calibration curve of temperature versus millivolt 
output. 


Description of Radiation Measurement 
and Data Reduction 

A summary of the radiation measurements made at 
lysimeter Y102C is given in table 7. All the sensors 
were purchased commercially, except for the 
miniature net radiometers, which were fabricated at 
the Water Conservation Laboratory at Phoenix, 
Ariz. The manufacturer's calibration had to be 
accepted for each instrument since no source of 
calibration was available. However, by comparing 
the output from similar sensois exposed during a 
period of steady solar radiation, a field calibration 
can be achieved. 

All the radiation sensors were located at sites that 
were mostly free from any obstructions above the 
plane of the sensing element. Fritschen net 
radiometers Nos, 432 and 567 and the inverted 
Eppley No. 7583DI were all mounted on the end of a 
2-meter boom that was pivoted from a 4-inch 
diameter mast 4 meters high. The mast was located on 
the north side of the lysimeter with the boom con- 
taining the sensors extending to the south over the 
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Table 6,-A summary of temperature and wind measurements made at Coshocton, Ohio, summer 1969 



Atr temper- 
ature (shield- 
ed and ven- 
tilated). 

Thermopile 

Instantaneous 
mV output at 

5 -min inter- 
vals converted 
to temperature 
(av for hour) 

L & N speedo- 
ina\ G, mV 
strip chart 
recorder. 

0.005^ C 

5 levels above 
crop surface. 
Sonic specific 
sensor levels 
given on daily 
profile data of 

(C).(G/ni3), 

(cm/s). 

Lysimeter Y102C. 

Wet-bulb 
temperature 
(shielded and 

- - do - - 

- - do - - 

- - do - - 

0.01 0^^ C 

- - do - - 
( 1 0 levels 
in all). 

- - Do. -- 

Watershed WS 101 
site center. 

ventilated) 

Air temt^er- 
aturc (shield- 
ed and venti- 
lated) 

- - do - - 

Instantaneous 
temperature at 
end of each 10- 
min period 
(av. for hour). 

L&N 

spcedoinax 

G tempera- 
ture strip 
chart recorder. 

0 25^ F 

5 levels above 
crop surface. 

Lysimeter Y 102C. 

Wet-bulb 
temperature 
(shielded and 

Thermo- 
couple with 
ceramic wick. 

-- do -- 

--do - - 

0.25** F 

- - do -- 

-Do.-- 

ventilated). 

Soil 

temperature. 

Thermo- 

couple. 

- - do - - 

- - do ~ - 

0.25** F 

0.05 m below 
ground surface. 

--Do.-- 

Soil 

surface 

temperature. 

Barnes 
infrared 
thermometer, 
IT-3, #493. 

Integrated 
hourly av. 
temperature 

Roy.son 

electronic 

integrator. 

o.r C 

1 ,0 m above 
soil surface. 

--Do.-- 

Wind 

1 Wong 

Av wind 

Sodcco 

1 count for 

2.4. 2 6, 3.0. 

Watershed WS 101 


cup 

anemometer 

speed for 
each hour. 

printing 

counter. 

e.ic]i 0.1- 
mi wind. 

and 2.0 111 
above crop 
surface. 

(all runs). 


2 Thorn- 

thwaite 

cup 

anemometers. 

-- do -- 

Mechanical 

counters 

(visual 

reading). 

1 count for 
each 1.5-111 
wmd 

2.4, 2.6, 3.0. 

3.8, 4.6, 5.4 m 
h - 2.4 m. 

Watershed 

WS 101 
(7/24/69). 


3 Beckman 
and Will tt ley 
cup 

anemometers. 

- - do - - 

Mechanical 

counters 

(microfilm 

recording 

15 min. 

- - do - - 

1.8. 2.3, 3.3 ni 
(10/14) 
h = 1.3 m. 

0.5, l.l, 1.6 m 
(10/16) 
h = 0 cm. 

Watershed 

WS 101 
(10/14), 

(10/16). 

Wmd 

direction. 

Beckman 
and Whittley 
vane. 

Compass 
point av. 
for hour. 

Strip chart 
recorder. 

±10° 

10 in above 
soil. 

Lysimeter Y102C. 



Table l.—A summary of radiation measurements made at Coshocton, OhiOt summer 


Observation | 

Sensor 

and 

senal no. 

Calibration 

constant 

Value 

recorded 

Recording 

system 

Sensitivity 
of recording 
system 

Height 
of sensor 

! 

Location 




^fv/langl€y 



Langley/mifi 

M 


Net 

radiation. 

1 Frit'jchen 
miniature 
net #432. 

3.12 

Integrated 
hourly from 
continuous 
measurement. 

Royson 

electronic 

integrator. 

0.005 

\ above 

crop 

surface. 

Lysimeter 

Y102C. 


2 Fntschen 
miniature 
ner ^567. 

3.42 

- - do - - 

- -do -- 

.0005 

- -do -- 

- - do • - 


3 Fntschen 
miniature 
net #548. 

3.43 

- - do ' - 

- - do - - 

.0005 

-- do -- 

Center of 
treated water- 
shed WS 101, 

Incoming 
total hemi- 
spherical 
radiation. 

Beckman &. 
Whitley 
total ven- 
tilated. 

1.27 

--do-- 

- - do - - 

.001 

- - do - - 

Lysimeter 

Y102C. 

Sod heat 
flux. 

Beckman & 
Whitley 
soil heat 
flux plates 
#540, #541. 

7.7 

--do-- 

- - do - - 

.00015 

0,05 below 
soil surface. 

--do-- 

Incoming 

solar 

radiation. 

Eppley 180" 
temperature 
compensated 
#7734. 

5.70 

Continuous 
trace and 
integrated 
hourly. 

L & N speedo- 
max H and 
integrator. 

,0015 

4 above 

crop 

surface. 

--do-- 

Reflected 

solar 

radiation. 

Eppley 
precision 
model 2 
#758301 
(inverted). 

5.35 

- - do - - 

L & N speedo- 
max H and 
disc inte- 
grator. 

.0015 

1 above 

crop 

surface. 

- - do - - 


lysimeter. Thus, the sensors were completely free of 
obstructions except for the mast, which, for practical 
purposes being to the north of the sensors, had no 
influence. 

Many radiation sensors were connected directly to 
electronic integrators. Each integrator system 
consisted of a Royson lectro-count which converted 
the incoming millivolt signal to a pulse rate, thepulse 
rate being directly proportional to the magnitude of 
the input signal. This pulse output from the electronic 
integrator was then totalized on Sodeco printing 
counters. The Sodeco units were programed to print 
the total counts for each hour on a paper tape and to 


reset to zero after printout. Thus, the integrated value 
of radiation from a sensor could be obtained by 
multiplying the count output for each hour by the 
integrator calibration constant and then by the sensor 
calibration constant. 

Two Leeds and Northrup Speedomax H strip chart 
recorders were equipped with disc integrators and 
General Electric printing demand meters. These units 
had been assembled by Eppley laboratories for use 
with their Eppley pyranometers. The output from the 
printing demand meter gave the total value of 
integrated radiation for each half-hour period in 
langleys. 
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APPENDIX C-AVERAGE HOURLY PROFILE DATA 


This appendix lists average hourly profile data for energy balance components at lysimeter and 
each day on which observations were made. Data are midwatershed sites. 

presented in the following order: 2. Smoothed profile data from log-linear plots of 

temperature, humidity, and wind at lysimeter and 
1, Average hourly lysimeter evaporation and midwatershed sites. 
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10 CM 0 = O.G CM CROP HEIGHT 



SMCOTHEO PROFILE DATA FRCM LOG-LINEAR PLOTS OF TE N'PERATURE t HUMIDITY, AND WIND 



o 

t 

% 

« 

• 

■ 

• 

* 

• 

• 

• 

• 

* 

• 

• 

tt 

• 

• 

• 

o 

o 


rg 

ITi 


nO 

o 

r- 

CM 




m 

m 

CO 

r- 

CM ro 

CO 

LU 

r-4 

tr\ 


C\J 

O' 

un 

'.n 

r\j 

O 

CM 

\X\ 


o 

r»H 


m 

O 

'4' 

c/0 


f«-4 







CM 

CM 

CM CM 


CM CM 

C\J 

CM 

CM 

CM 

■V 




















5: 




















o 




















o 





















o 

• 

• 

• 

« 

* 


• 

• 

f 

• 

• 

• 

• 

• 

• 


• 

• 


-H 

o 


fsJ 

“O 

CO 

vjD 

iTN 


sO 

o 


o 

M3 

O' 

o 

pm«4 

O' 

Ch 

X 


r- 

o 

-C 

't 

r\t 

h- 

vO 

o 

rA 

fsj 


o 

O 

o 

CM 

o 

o 

O 











1— 1 

rA 





rA 





























ro 






















X o 

00 

vO 

fO 

O 

o 

lA 

in 

CM 

00 

o 

CM 


lA 

r- 

lA 

A 

LA 

A 




^ o 

(r\ 

cr 

tA 

o 

rH 

ro 

xO 

O' 

O' 

O' 

QC 

O 


h- 

A 

O 

A 




X rA 

• 

• 

• 

• 

• 

« 

» 

• 

• 

• 

• 

« 

• 

« 

• 

9 

9 

9 




e> 

7 

7 


CO O' 

c^ 

o 

O' 

CT' 

6 

6 

O' 

O 

o 

O' 

00 

CO 



a 


*«>• 













pH 

pH 





9 

bxri 





















S' 

X 


> 



















u 

o 


>- 






















K»-i 



















o 



C 



















o 

o 


HH O 

<r 

u 

<A 

rn cj 

A 

m 

o 

CM 


CM 

A 

A 

4" 

u> 

<r 

4* 

O 

rH 

h— 


SI *A 

lA 

<r 

O 

4^ 

o 

O 

in 

cc 

lA 

CM O 

r— J 

ro 

A 

CM 

A 

O' 

h- 


u 


:j 

■ 

t 

• 

• 

« 

• 

• 

• 

• 

• 

• 

• 

« 

« 

• 

9 

9 

9 

O 

a 


X 

r- 

r- 

CO 

O' 

o 

o 

o 

O' 

O 

pH 

CM 

(M 

fA 

pH 

o 

O' 

CC 

00 


X 







♦H 

1-H 



pH 

pH 

pH 

pH 

pH 

pH 

pH 




< 



« 




















O 


LO 



















U) 

o 


CU 



















pH 



< 






















o 

O' 

O' 


r- o 


r- 

O' 

00 

o o 

ro 

lA 

A 

CC 

C\J 

4' 

CM 

II 



o 

00 


LA 

O' 

4^ 

0> 

fA 

(A 

pH 

A 00 

A 

CM 4- 

o 

LA 

00 

cr 





• 

• 

• 

• 

• 

« 

• 

• 

♦ 

• 

• 

• 

» 

9 

9 

9 

9 

9 

00 

o 



a; 

00 

CT 

o 

CM 

rH 

1^ 


CM 

m 

<r 

tA 

lA 

4- 

rn 


o 

o 

A 

1 







rA 

rH 

1-^ 

pH 

pH 

1— • 

pH 

pH 

pH 

r-< 

pH 

P-H 

A 

pH 

2 

CM 

O' 




















UJ 

o 

>o 




















X 

rA 


UJ 



















uu 

V 

N*. 

cc 



















cc 


CM 

X 



















X 

UJ 

CM 

1- o 

vO 

LA 

CM 

rg 

o 

r- 

4* 

O' 

GO 

(A 

CM 

1— J 

in 

00 

A 

cn 

4^ 

A 

»- 


< r-l 

O 

o 

n 

00 

CO 

in 

rA 

'Tk 

lA 

o 

4* 

O' 

O' 

o 

A 

A 

00 

4' 

< 


V 

cc 


« 


• 

• 

• 

• 

• 

f 


• 

• 

• 

t 

9 

9 

9 

9 

LU 

oo 

lA 

UJ 

CD 

00 

CA 

rJ 

CM 

CM 

»— 4 

p-H 

CM 

CA 

u^ 

lA 

4^ 

4- 

CM 

<p-I 

o 

o 

X 



Q. 






tA 

rA 

rH 


i-H 

A 

.H 

pH 

f-H 


pH 

pH 

pH 


iX 


3C 



















LL 

UJ 

1 

UJ 



















o 



J- 




















UJ 

LU 




















H 


H 




















X 


<r 

o 



















O 

i/> 

o 

Q C 



















A 

> 


2: »- 



















UJ 

— 1 


U) oc 

in 

>0 


00 

O' 

o 

rA 

CM 

fO 

4* 

m 

A 

h" 

CO 

a 

o 

pH 

CM 

X 



Uj 

UJ Cl- 







tA 

pH 



i-i 

p-H 

i-H 

fH 


CM 

A4 

AJ 




a. 






















M ll 





* 

















K O 





















40 
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SNCCTHED profile DATA FRCM LQG-LINEAR PLOTS OF TEVPERATURE, HUMIDITY, AND WIND 
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SMCOTHEO PPOFILE DATA FRCM LOG~LINEAR PLOTS OF TEMPERATURE, HUMIDITY, AND WIND 
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SMCGTHED PROFILE DATA FROM LOG-LINEAR PLOTS OF TEMPERATURE, HUMIDITY, ANC V,IND 
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SMOOTHED PROFILE DATA FROM LOG-LINEAR PLOTS OF TEMPERATUREt HUMIDITY» AND WIND 
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SMOOTHED PROFILE DATA FROM LOG-LINEAR ^.,:TS OF TEMPERATURE, HUMIOETY, AND WIND 
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SMCOTHED PROFILE DATA FROM LOG-LINEAR PLOTS OF TEKPERATUREf HUMIDITYt AND WIND 
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APPENDIX D---AVERAGE HOURLY EVAPORATION 


This appendix reports average houily evapotrans- 
piration by the various meteoiological equations 
using the smoothed profile data and average houily 
radiation data. Also presented aie the computed 
Bowen Ratio (B) and Richardson’s Number (Rl). 


The values in the probable erim columns gave the 
range of piobuble enor that can be expected in each 
computed value on the basis ot instrument accuracy 
and sensitivity. Missing data arc due to equipment 
malfunctiom as indicated by an asterisk. 
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